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1.  Introduction 

Digitized  terrain  data  corresponding  to  a  given  geographical  area  are  generally 
available  in  a  rectangular  grid  at  some  fixed  resolution  (i.e.,  distance  between  neighboring 
(x,y)-coordinates).  In  the  modeling  and  simulation  of  processes  that  interact  with  terrain, 
it  is  often  true  that  the  higher  the  resolution,  the  better.  Thus,  it  is  often  desirable  to  be 
able  to  construct  a  model  of  a  piece  of  terrain  that  has  realistic  terrain  features  that  are 
smaller  than  the  resolution  at  which  the  digitized  data  are  available.  We  must  use  the 
term  "realistic"  rather  than  "real"  since  it  is  impossible  to  resolve  any  actual  terrain  fea¬ 
ture  smaller  than  the  resolution  of  the  data.  For  the  purposes  of  this  report,  a  "micro¬ 
terrain  feature"  is  any  terrain  feature  that  is  smaller  than  that  fixed  resolution,  and  "micro- 
terrain"  is  a  collection  of  such  micro-terrain  features.  A  more  precise  mathematical  defi¬ 
nition  of  a  terrain  feature  and  the  measure  of  its  size  is  given  in  section  2. 

In  a  previous  report,  the  author  described  an  algorithm  that  creates  a  smooth  Vari¬ 
able  Resolution  Terrain  (VRT)  surface  that  matches  a  set  of  digitized  terrain  altitudes 
(DTAs)  to  within  a  prescribed  tolerance  at  those  (x,  y)-coordinates  at  which  the  DTAs  are 
defined.1  In  this  report,  that  methodology  is  applied  to  the  problem  of  adding  micro¬ 
terrain  to  a  set  of  DTAs. 

The  VRT  model  is  described  in  section  2.  The  process  of  fitting  a  VRT  surface  to  a 
set  of  DTAs  is  explained  in  section  3.  In  section  4,  the  basic  methodology  for  adding 
micro-terrain  is  explained,  with  an  example.  Section  5  contains  some  thoughts  on  model 
validation  procedures  and  future  work. 

2.  The  VRT  Model 

The  VRT  model,  developed  at  the  U.S.  Army  Ballistic  Research  Laboratory  (now 
the  U.S.  Army  Research  Laboratory)  in  the  early  1990s,  is  a  model  that  represents  basic 
topography  as  a  continuous  surface,  capable  of  being  viewed  at  any  desired  resolution.2 
This  surface  is  the  superposition  of  individual  terrain  features,  or  "hills"  for  simplicity, 
each  of  which  is  described  by  a  closed  form  mathematical  function  that  is  continuous 
everywhere  and  has  continuous  partial  derivatives  almost  everywhere.  While  examples  of 
many  different  types  of  terrain  have  been  constructed  using  the  VRT  model,  they  are  all 
"generic"  in  the  sense  that  they  were  developed  from  basic  principles  and  do  not  represent 

1  Wald,  Joseph  K.  ", Solving  the  Inverse  Problem  in  Terrain  Modeling ."  ARL-TR-605, 

U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  October  1994. 

2  Wald,  Joseph  K.  and  Patterson,  Carolyn  J.  "A  Variable  Resolution  Terrain  Model  for 
Combat  Simulation."  BRL-TR-3374,  U.S.  Amy  Ballistic  Research  Laboratory, 

Aberdeen  Proving  Ground,  MD,  July  1992. 
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any  particular  piece  of  Earth  topography.  For  generic  terrain  creation,  the  distribution  of 
the  sizes  of  the  hills  is  based  on  the  idea  of  self-similarity  (i.e.,  invariance  with  respect  to 
scale),  which  is  embodied  in  the  power  law 

D  =  Ks~\  (2.1) 

in  which  j  is  a  dimensionless  scale  factor  associated  with  a  hill,  AT  is  a  constant  that 
depends  on  terrain  type,  and  D  is  the  areal  density  of  hills  as  a  function  of  s.  The  choice 
of  exponent  in  this  power  law  ensures  self-similarity  in  the  density  of  terrain  features. 
Integrating  this  power  law  produces  the  cumulative  distribution  function: 

t  «  s 

J  Kt~2dt  =  K(1  /  -  1/s).  (2.2) 

f  *nrin 

There  is,  in  practice,  a  range  of  scales,  smaJ,  for  which  the  power  law  holds.  Thus, 
to  build  each  hill,  it  is  only  necessary  to  draw  a  uniform  random  number  u  from  the  inter¬ 
val  (0, 1  —  [jmjn  /  ^max]),  with  the  hill  scale  factor  being  s  =  sm;n  /  (1  —  u).  The  location  of 
the  hill  is  randomly  chosen  in  the  desired  area.  The  complete  terrain  surface  is  defined  by 
the  superposition  of  all  of  the  hills,  i.e., 


N 


T(x,y)=^fk(x,y). 


k= 1 


The  form  of  a  single  hill  function,  fk(x,  y),  is  given  by: 


fk(x,  y)  =  sk  hkexp(  - 


“T  [<**£*(*, y)  +  (l  -ak)Mk(x,y) 
skPk  L  J 


°k 


where 


ax  =f*-^*-^-Jcos2  Xk, 
=  sin  2A*’ 


(2.3) 


X  (2.4a) 


r  - 

Ek(x>  y )  =  J^iC*  - 1*)2  -  a2(x  -  |*)(y  -  Tjk)  +  a3(y  -  ;7*)2J2,  (2.4b) 


(2.4c) 


(2.4d) 


and 


2 


and 


<*3  =  £k-  ^  -  j-1  sin2  A*, 


(2.4c) 


Mk(x,  y )  =  max  (lAj  (x,  y)  |, \A2(x,  y)|),  (2.4f) 

where 

y)  =  £*(*  -  lit)  cos  A*  +  —  (y  -  77*)  sin  A*,  (2.4g) 

and 


>)  =  -(>’- Tjk) cos  A*  -£*(*-  £*)  sin  A*,  (2.4h) 

This  formulation  is  a  substantial  revision  and  extension  of  the  original  hill  function 
definition.3  Varying  the  parameters  ak,  |*,  rjk,  hk,  pk,  ek,  Xk,  ak,  and  sk  produces  hills  in 
a  variety  of  sizes  and  shapes.  While  all  of  the  parameters  interact  to  some  extent,  roughly 
speaking,  the  parameters  |*  and  rjk  govern  the  location  of  the  "center"  of  the  hill,  hk,  pk, 
and  sk  govern  the  "height"  and  "width"  of  the  hill,  ek  governs  the  "eccentricity"  of  the 
hill.  A*  governs  the  "orientation"  of  the  hill,  and  ak  governs  the  "slope"  of  the  hill.  The 
parameter  ak ,  whose  value  always  lies  in  the  interval  [0, 1],  governs  the  shape  of  the  hori¬ 
zontal  cross  section  of  the  hill.  When  ak  =  1,  this  cross  section  is  elliptical.  When 
ak  =  0,  the  cross  section  is  rectangular.  For  other  values  of  ak,  the  cross  section  has  a 
more  complex  geometry.  Rectangular  hills  are  useful  in  the  modeling  of  houses,  roads, 
and  other  man-made  objects  as  integral  elements  of  the  terrain. 

3.  Constructing  a  VRT  Surface  that  Fits  a  Set  of  DTAs 

The  process  of  fitting  a  VRT  surface  to  a  set  of  DTAs  has  been  documented.4  How¬ 
ever,  since  that  methodology  plays  an  important  part  in  the  process  of  adding  micro¬ 
terrain,  it  is  appropriate  to  briefly  describe  it  here. 

The  problem  is  to  construct,  in  closed  form,  a  surface  z  =  F(x,  y),  with  the  proper¬ 
ties  that  F(x,  y)  is  continuous  everywhere,  has  continuous  partial  derivatives  Fx(x,  y)  and 
Fy(x,  y)  almost  everywhere  (i.e.,  everywhere  except  a  set  of  measure  zero),  and  that 
"matches"  a  given  set  of  DTAs  to  within  a  prescribed  tolerance,  S  (i.e.,  for  each  DTA 

3op.  tit. 

4  Wald,  Joseph  K.  ", Solving  the  Inverse  Problem  in  Terrain  Modeling ."  ARL-TR-605, 

U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  October  1994. 


3 


point  ( x,y,d ),  \F(x,y)-d\  <,8).  [The  advantage  of  obtaining  such  an  F(x,y)  is  clear. 
At  virtually  any  point  in  the  geographical  area  under  consideration,  both  the  terrain  alti¬ 
tude  and  the  slope  of  the  terrain  in  any  direction  are  obtained  simply  by  evaluating 
F(x,  y)>  F x(x>  y)y  and  Fy(x,  y)  at  the  desired  (x,  y)-coordinate.] 

First,  the  DTAs  are  detrended  by  fitting  a  plane  through  them  and  subtracting  the 
height  of  this  plane  from  each  of  the  DTAs.  This  produces  a  set  of  "residual"  DTAs. 
Next,  the  residual  DTA  (x,  y,  d)  of  largest  magnitude  is  found  and  a  hill  is  constructed, 
the  center  of  which  is  at  or  near  ( x ,  y)  and  the  height  of  which  is  equal  to  or  close  to  d. 
An  iterative  loop  is  used  to  compute  the  values  of  the  VRT  parameters  in  the  construction 
of  this  hill.  Here,  the  best  hill  is  the  one  that  leaves  the  smallest  modified  residual  DTAs 
in  the  neighborhood  of  ( x ,  y)  when  the  hill  is  subtracted  from  the  residual  DTAs.  This 
hill  is  subtracted  from  each  of  the  residual  DTAs  to  create  a  new  set  of  residual  DTAs. 
The  new  largest  residual  DTA  is  located  and  the  hill  fitting  process  is  repeated  until  all  of 
the  residual  DTAs  are  smaller  in  absolute  value  than  8.  The  VRT  "fitting"  surface, 
F(*>  y)>  is  then  defined  to  be  the  sum  of  the  fitting  plane  and  all  of  the  fitting  hills, 

Note  that  since  a  pseudo-random  number  generator  is  used  in  the  iterative  loop, 
simply  changing  the  random  number  seed  will  produce  a  slight  variant  for  F(x,y).  In 
fact,  there  are  an  infinite  number  of  surfaces  that  satisfy  the  conditions  for  F(x ,  y). 

4.  Procedure  for  Adding  Micro-Terrain 

When  constructing  a  generic  VRT  surface,  micro-terrain  is  included  simply  by  set¬ 
ting  sm j,,  equal  to  the  scale  of  the  smallest  desired  terrain  feature  (in  equation  2.2).  How¬ 
ever,  adding  micro-terrain  to  a  set  of  DTAs  is  a  more  complicated  matter.  The  first  step  is 
to  construct  a  VRT  surface,  Fj(x,y),  that  fits  the  DTAs  to  within  a  prescribed  tolerance, 
8.  Next,  empirical  cumulative  distributions  are  formed  for  the  VRT  parameters  ak,  hk, 
Pk,  ok,  and  sk  from  the  parameters  of  the  hills  used  to  construct  F,(x,y).  The 

proportionality  constant  K  (from  equations  2.1  and  2.2)  is  estimated  from  the  cumulative 
distribution  of  the  values  of  the  sk  parameter  by  fitting  a  curve  of  the  same  form  as  equa¬ 
tion  2.2  to  these  values.  Here,  the  value  for  is  just  the  scale  of  the  smallest  terrain 
feature  used  to  construct  Fi  (x,  y). 

Having  approximated  the  elements  of  the  VRT  model,  the  next  step  is  to  specify  a 
new  denoted  by  s*^,  corresponding  to  the  scale  of  the  smallest  micro-terrain  fea¬ 
ture  desired,  and  then  create  a  set  of  micro-hills  by  successively  drawing  uniform  random 
numbers  «  from  the  interval  (0, 1  -  [s* ^  /  j^]),  with  the  scale  factor  of  a  given  hill 
being  /  (1  -  «).  The  location  of  each  hill  is  randomly  chosen  in  the  area  covered  by 
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the  VRT  surface.  The  complete  micro-terrain  surface  is  defined  by  the  superposition  of 
all  of  the  micro-hills,  i.e., 


at, 

G(x,  y)  =  2  mk(x>  y)>  (4.1) 

*=i 

where  G(x,y )  is  the  "micro-surface"  and  each  of  the  micro-hills,  mk(x,y),  has  the  same 
form  as  fk(x,  y)  in  equation  2.4,  with  the  parameters  ak,  hk,  pk,  ek,  Xk,  and  ak  randomly 
selected  from  the  respective  empirical  cumulative  distributions. 

At  this  point,  the  sum  F,(jc,  y)  +  G(x,  y)  is  not  formed.  While  it  is  true  that  this 
sum  does  approximate  the  DTAs  (to  some  degree)  and  also  contains  micro-terrain,  the 
condition  |Fi(x, y)  +  G(x, y)-d\<S  may  be  violated  for  some  (*, y, d)  (i.e.,  the  fit  of 
F\(x, y)  to  the  DTAs  may  have  been  disturbed  by  the  addition  of  the  micro-terrain). 
Instead,  Ft(x,y)  is  discarded  and  a  second  fitting  surface  is  constructed,  although  in  a 
slightly  different  manner.  Instead  of  the  flat  plane,  z  =  0,  construction  starts  with  the 
micro-surface,  z  =  G(x,y),  and  to  this  surface  a  "detrending"  plane  and  a  sequence  of 
empirically  determined  hills  is  added,  using  the  same  type  of  iterative  loop  as  described 
in  section  3.  The  product  of  this  loop  is  a  new  approximating  surface,  F2(x,  y).  Since  the 
loop  termination  condition  is  the  same  as  before,  the  surface  S(x,  y)  =  F2(x,  y)  +  G(x,  y) 
satisfies  the  conditions  for  a  fitting  surface  and  also  contains  micro-terrain  down  to  the 
desired  scale. 

This  procedure  works  well  when  the  DTAs  represent  a  single  terrain  type.  When 
the  DTAs  are  extensive  enough  to  include  more  than  one  terrain  type,  however,  the  empir¬ 
ical  cumulative  distributions  should  really  be  localized  (i.e.,  the  selection  of  parameters 
for  the  micro-hills  should  be  based  on  the  macro-hills  in  the  same  area).  Also,  the  pro¬ 
portionality  constant  K  may  need  to  be  localized. 

Figures  1  through  3  illustrate  the  procedure  for  adding  micro-terrain  to  a  set  of 
DTAs.  A  set  of  DTAs  representing  a  25-meter  square  of  terrain  from  the  National  Train¬ 
ing  Center,  Fort  Irwin,  California,  was  selected  from  a  larger  set  of  DTAs  provided  by  the 
U.S.  Army  Topographic  Engineering  Center.5  The  grid  spacing  for  the  DTAs  is  5  meters, 
resulting  in  a  total  of  36  DTAs  in  the  set  Figure  1  depicts  the  grid  of  DTAs.  Clearly,  any 
actual  terrain  feature  contained  in  the  piece  of  terrain  from  which  the  DTAs  were 
extracted,  but  smaller  than  the  grid  size,  cannot  be  "captured"  by  the  DTAs.  In  figure  2,  a 
VRT  surface,  constructed  using  the  methodology  described  above,  has  been  superim¬ 
posed  on  the  DTAs.  The  maximum  deviation  of  the  VRT  surface  from  the  DTAs  (at 

*F.  Raye  Norvelle,  U.S.  Army  Topographic  Engineering  Center.  Private  Communication, 

March  1994. 
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Figure  2.  A  25-meter  square  -  VRT  surface. 


Figure  3.  Close-up  of  a  5-meter  square. 
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those  (x,y)-coordinates  at  which  DTAs  were  collected)  is  8  centimeters  and  the  average 
deviation  is  4  centimeters.  A  total  of  939  hills  were  constructed  in  the  combined 
"micro/fitting"  process,  with  micro-hills  being  constructed  down  to  a  scale  of  1  centime¬ 
ter.  The  construction  of  the  VRT  surface  took  slightly  less  than  2  seconds  on  a  CRAY 
YMP  computer.  Figure  3  shows  some  of  the  smaller  micro-hills  in  an  "ant’s-eye"  view  of 
a  5-meter  square  subset  of  the  VRT  surface. 

5.  Model  Validation  and  Future  Work 

The  fact  that  the  VRT  surface  matches  the  DTAs  to  within  the  prescribed  tolerance 
is  easily  checked  by  direct  calculation.  Since  S(x,  y)  is  a  finite  sum  of  continuous  func¬ 
tions  with  continuous  derivatives  of  all  orders  almost  everywhere,  it  inherits  these  proper¬ 
ties. 

However,  the  claim  that  the  micro-terrain  features  are  realistic  is  harder  to  verify, 
especially  since  the  term  "realistic"  is  not  defined  mathematically.  The  procedure  for 
adding  micro-terrain  clearly  depends  strongly  on  the  self-similarity  assumption  underly¬ 
ing  the  basic  theory  for  VRT.  This  procedure  works  as  long  as  the  hills  constructed  lie 
within  the  range  of  scales  for  which  that  assumption  is  valid.  Unfortunately,  it  is  impossi¬ 
ble  to  know  in  advance  what  that  range  is.  There  is,  however,  an  important  clue  in  the 
procedure  that  can  help  shed  some  light  on  this  matter.  Recall  that  the  proportionality 
constant  K  is  estimated  from  the  empirical  cumulative  distribution  of  the  values  of  the  sk 
parameter.  The  degree  to  which  that  empirical  distribution  matches  the  theoretical  distri¬ 
bution  of  equation  2.2  gives  a  measure  of  how  closely  the  real  piece  of  terrain  from  which 
the  DTAs  were  extracted  matches  an  "ideal"  VRT  surface.  If  that  match  is  good,  then  the 
assumption  of  self-similarity  is  valid,  at  least  down  to  the  scale  corresponding  to  the  reso¬ 
lution  of  the  DTAs.  If  the  match  is  not  very  good,  there  may  still  be  a  restricted  range  of 
scales  for  which  a  good  match  can  be  found. 

In  any  case,  to  validate  the  micro-terrain  creation  model,  it  is  necessary  to  define  a 
process  that  demonstrates  that  for  a  given  set  of  DTAs,  the  micro-terrain  reflects,  in  some 
way,  the  small  terrain  features  resident  in  the  piece  of  terrain  from  which  the  DTAs  were 
extracted.  Of  course,  since  DTAs  are  isolated  points,  it  is  impossible  to  malcp.  a  direct 
comparison.  In  addition,  it  must  be  remembered  that  the  model  does  not  claim  that  it  is 
recreating  actual  small  terrain  features  (which  are  unknown  below  the  resolution  of  the 
DTAs),  but  rather  that  the  micro-terrain  it  is  creating  is  representative  of  that  type  of  ter¬ 
rain.  These  facts  suggest  that  some  type  of  statistical  comparison  is  in  order.  For  exam¬ 
ple,  for  a  set  of  DTAs  collected  at  a  given  resolution,  it  is  possible  to  construct  a  VRT 
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surface  with  micro-terrain,  S(x ,  y),  using  a  coarser  subset  of  the  DTAs.  In  fact,  it  is  pos¬ 
sible  to  construct  many  such  surfaces,  S,(x,  y),  i  =  1, 2, . . . ,  L.  Then,  for  each  DTA  point, 
(x,y,d),  not  used  in  the  fitting  process,  compute  the  empirical  distribution 
—  St(x, y), i  =  1,2 ,L.  This  provides  a  basis  for  comparison  between  the  stochastic 
model  and  the  data  at  the  points  for  which  a  comparison  is  possible.  At  this  point  there 
are  many  different  approaches  that  could  be  used  to  make  the  comparison.  These  will  be 
investigated  in  future  work.  It  should  be  pointed  out  that  this  validation  methodology  can 
also  be  considered  as  part  of  a  design  feedback  loop,  in  which  the  characteristics  of  added 
micro-terrain  can  be  adjusted  to  match  a  given  type  of  terrain.  Although  this  process  is 
still  limited  to  the  resolution  of  the  DTAs,  it  can  some  give  some  measure  of  confidence 
when  extrapolating  to  higher  resolution. 

By  repeating  the  micro-terrain  creation/validation  process  for  large  number  of  sets 
of  DTAs  from  a  given  geographical  area,  it  should  eventually  be  possible  to  characterize 
terrain  types  in  terms  of  the  parameters  of  the  VRT  model.  The  undertaking  of  such  a 
massive  task  is,  however,  well  beyond  the  scope  of  the  present  work. 


8 


NO.  OF 

COPIES  ORGANIZATION 


2  ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CTR 
ATTN  DTIC  DDA 
CAMERON  STATION 
ALEXANDRIA  VA  22304-6145 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TP 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 


ABERDEEN  PROVING  GROUND 

5  DIR  USARL 

ATTN  AMSRL  OP  AP  L  (305) 


NO.  OF 

COPIES  ORGANIZATION 


1  HQDA  SARD  TR  MS  K  KOMINOS 
WASH  DC  20310-0103 

1  HQDA  SARD  TR  DR  R  CHAIT 
WASH  DC  20310-0103 

1  US  ARMY  CORPS  OF  ENGINEERS 
WATERWAYS  EXPERIMENT  STATION 
MOBILITY  SYSTEMS  DIVISION 

ATTN  CEWES  GM  ANDREW  W  HARRELL 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  39180-0631 

2  US  ARMY  TOPOGRAPHIC  ENG  CTR 
ATTN  CETEC  PD  DA 

RAY  NORVELLE 
GEORGE  LUKES 
7701  TELEGRAPHIC  RD 
ALEXANDRIA  VA  22310-3864 

1  INSTITUTE  FOR  DEFENSE  ANALYSIS 
ATTN  WARREN  K  OLSON 
ROOM  242W 

1801  N  BEAUREGARD  ST 
ALEXANDRIA  VA  22311-1772 

ABERDEEN  PROVING  GROUND 

108  DIR  USARL 

ATTN  AMSRLWTW  CHARLES  H  MURPHY 
AMSRL  WT  WE 
JUDITH  TEMPERLEY 
MARY  ANNE  FIELDS 
CAROLYN  PATTERSON 
JOSEPH  WALD  (100  CP) 

AMSRL  Cl  SB 
VIRGINIA  KASTE 
THOMAS  KENDALL 
CHARLES  HANSEN 
TERRY  PURNELL 

2  DIR  USAMSAA 

ATTN  AMXSY  ED  JACK  MEREDITH 

AMXSY  CD  DWAYNE  NUZMAN 


10 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-866 _ Date  of  Report  September  1995 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate.  _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ _ _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail ) 
(DO  NOT  STAPLE) 


